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ABSTRACT 

The state of the hot gas in clusters of galaxies is investigated with a set of 
model clusters, created by assuming a polytropic equation of state (P = 1.2) 
and hydrostatic equilibrium inside gravitational potential wells drawn from a 
dark matter simulation. Star formation, energy input, and nonthermal pressure 
support are included. To match the gas fractions seen in non-radiative hydrody- 
namical simulations, roughly 5% of the binding energy of the dark matter must 
be transferred to the gas during cluster formation; the presence of nonthermal 
pressure support increases this value. In order to match X-ray observations, 
scale-free behavior must be broken. This can be due to either variation of the 
efficiency of star formation with cluster mass M500, or the input of additional 
energy proportional to the formed stellar mass Mp. These two processes have 
similar effects on X-ray scalings. If 9% of the gas is converted into stars, inde- 
pendent of cluster mass, then feedback energy input of 1.2 x 10~^Mpc'^ (or ~ 1.0 
keV per particle) is required to match observed clusters. Alternatively, if the 
stellar mass fraction varies as oc M^olJ-^^ then a lower feedback of 4 x lO'^Mpc'^ 
is needed, and if the stellar fraction varies as steeply as oc Mg^o-^^ then no addi- 
tional feedback is necessary. The model clusters reproduce the observed trends 
of gas temperature and gas mass fraction with cluster mass, as well as observed 
entropy and pressure profiles; thus they provide a calibrated basis with which to 
interpret upcoming SZ surveys. One consequence of the increased gas energy is 
that the baryon fraction inside the virial radius is < 90% of the cosmic mean, 
even for the most massive clusters. 

Subject headings: cosmology: theory — galaxies:clusters:general — intergalactic 
medium — X-rays:galaxies:clusters 
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Introduction 



In the near future, clusters of galaxies will prove to be an increasingly important probe 
into the abundance and properties of both dark matter and dark energy. Data on large 
numbers of clusters are becoming avail able from surveys at many wavelengths: in the opti- 



cal/infrared, from galaxy overdensities (iKoester et al.ll2007l : lYee et al.ll2007l : lEisenhardt et al. 



2008 : iLopes et al.l 120091) or from weak gravi t ational lensing of background objects (IDahle 



20071): in the X-r av (IBohringer et al.l 120071 : iBurenin et al.l 120071 : iFinoguenov et al. 



2007 



Pacaud et al.ll2007l ) ; and in the microwave, by locating temperature decrements caused by the 



Sunyaev-Zeldovich (SZ) effect (IDobbs et a. 



20071 : iMuchovei et al.ll2007l : iLin et al. 



2008 



2006 



Kosowskv 



2006 



Bonamente et al 



Malte Schafer &: Bartelmann 



2008 



Staniszewski et al.ll2008l ). 



Properly interpreting such surveys will require a sound understanding of the state of 
the Intra-Cluster Medium (ICM). Based on the assumption that the thermal energy in the 
gas comes solely from gravitational collapse, one may deri ve expected self-similar scalings 
between halo niass a nd temperature. X-ray luminosity, etc. (lKaiserlll99ll : iBryan &: Norman 
19981 : lEke et al.lll998r). However, thes e scahngs do not match the actual behavio r of observed 



clust e rs (e.g 



Arnaud fc Evrard 



2008; Sun et al 



2009 



1999: 



Vikhlinin et al. 2006; Arnaud et al. 2007; Pratt et al. 



Vikhlinin et al. 



20091 ) . Additional modification of the thermal state of 
the gas from non-g ravitational s ources will bring these predictions into better agreement with 
observed scalings (jKaiseiill99ll ). One such modification is star formation; by removing gas 
with the shortest cooling times, star formation leaves behind gas with a higher mean entropy 
( jVoit &: Bryarul200ll : iTozzi fc Normanll200ll : IVoit et al.ll2002l ) . Energy injection from preheat- 
ing or stellar and AGN feedback like wise affects scalings, as numerous simulations which in- 



clude such feedback haye shown fe.g. lEttori et al. 



200 



200 



Nagai et al. 



Dave et al. 



2007 



2006 



20081 : IPuchwein et al.ll2008l : 



Bqrgani et al 



2006 



Muanwong et al. 



Kav et al.ll2007l:ISiiacki et al.ll2007l: lBha.ttacharva et ahlbood : iBorgani et al 



Short fc Thomaj 20081 . and references therein) . 



The importance of these processes can be investigated analytically by co nsidering hydro- 



static, and usually i 


polytropic, gas residing in various c 


ark matter potem 


tials 


Suto et al. 


1998 




Balog 


1 et al. 


1999; 


Wu et allboool: Loewenstein 200ol: 


Tozzi & Norman 


2001; 


Komatsu & Se 


[jak 


2001 




Babul et al. 


2OO2I; 


Voit et al. 




2OO2I: 


Dos Santos & Dorel I2OO2I: . 


A.scas 


ibar et al. 


2003 




Shimizu et al 


2004 


; LaDi et al.ll2005 


Afshordi et a 


I2OO5I: Solanes et al.ll2005 


: Ascasibar & Diego 


2008 




Bower et al. 


20081; 


Ciotti & Pellegrini 


2008: I 


Vloodlev et al. 2008). 



Qstriker et al.l (120051 . hereafter Paper I) and iBode et al.l (l2007l . Paper II) developed a 
simple method for populating the halo Dark Matter (DM) potentials from N-body simula- 
tions with gas, based on the assumptions of hydrostatic equilibrium and a polytropic equation 
of state. This method can be calibrated against local X-ray observations, and then be used 
to create realistic sky maps (ISehgal et al.l 120071 ). Rather than assuming a simple analytic 
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profile, the distributions of lialo concentration, substructure, and triaxiality are included, as 
well as any spatial correlations, since the halos are drawn from an accurate large scale dark 
matter simulation. This procedure has advantages over performing full hydrodynamic sim- 
ulations. Because less computational expense is required, large volumes can be simulated. 
Also, one can vary parameters, such as the amount of star formation and feedback, without 
having to redo an entire simulation. 



This paper presents a number of improvements to the model presented in iPaper t and 



Paper 111 . A more sophisticated description of star formation is inclu ded, whic h accounts for 
gas recycling and allows for a calculation of metallicity evolution. iPaper III assumed that 
every cluster hal o has the same stellar mass fraction; here we allow this fraction to vary with 
halo mass. Also iPaper 111 did not allow for transfers of energy from the dark matter to the 
gas. It is evident from simulations that this does occur, so we allow for this transfer here. 



Because of these limitations, iPaper 111 required a large amount of feedback energy to match 
X-ray observations. As we shall see, the improvements presented here significantly reduce 
the feedback required. 

With this method, we examine how star formation and different forms of energy input 
affect the ICM. It is possible to produce model clusters which match existing observational 
constraints, and thus make predictions for upcoming surveys. The method is described in 
Sec. |2] and Appendix \M In Sec. [3] it is then applied to a simulation reflecting the current 
best-fit cosmological parameters and the results are compared to X-ray data; in particular 
the effects of varying star formation and feedback are examined. In Sec. IHa standard model 
is compared to further X-ray and SZ observation; we summarize and discuss implications in 
Sec. El 



2. Method 

To create a mock catalog, cluster-sized halos are extracted from a dark matter only 
N-body simulation, and then the gas distribution inside each halo potential is determined. 
This section describes the gas model and the construction of the halo catalog which is used. 



2.1. The Gas Prescription 

A given DM halo is enclosed in a cubic mesh, using a cell size, I, twice the N-body 
particle spline softening length. The DM density, pDfc, and gravitational potential, 0^, are 
computed for each cell k using a Particle-Mesh code with non-periodic boundary conditions. 
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The cell with the lowest potential, (po = MIN{(l)k), is considered the center of the cluster. 
Particle velocities in the rest frame of the cluster (defined as the mean velocity of the 125 
particles closest to the cluster center) are used to find t he KE per unit volume, ^tofc, in each 



cell. The virial radius (from spherical top-hat collapse; iBryan fc Normanlll998l ) can then be 



computed, as well as the halo virial mass, Md, and energy, Ed. 

It is assumed that gas originally followed the DM inside the virial radius, with density 
fcPDk and KE /c|tDfc, where /c = flb/flm- A certain amount of the gas, M^, will have turned 
into stars. We parameterize this mass by the fraction /* = M^/{fcM]:,) at z=0; the evolution 
of with redshift is described in Appendix |X1 This stellar mass presumably formed from 
the portion of the gas with the lowest entropy and shortest cooling time, so cells are ranked 
by binding energy and then are checked off, starting with the most bound, until the sum of 
their gas masses equals M^] the gas mass and energy in these cells are set to zero. Thus the 
initial mass Mg and energy Eg of the remaining gas are: 

Mg = Y,fcPTykl' , (1) 

k 

Eg = ^fc {0fepDfe + itofe} , (2) 

k 

summing over all cells inside r^r except those marked off for star formation. The surface 
pressure Ps on the gas at the virial radius is estimated from the kinetic energy in a buffer 
region eleven cells thick: 

. Nb 

summing over the N}, cells in the buffer region r^ij. < < r^r + 11/. 

From this initial state, it is assumed the gas is rearranges itself into hydrostatic equilib- 
rium in the DM potential, with a pol ytropic e quation of state of index T = 1.2 (as indicated 



by both observation and simulation; iPaper t ). The resulting gas pressure P and density p 



r 1 

are given by Pk = Po^k"^ and pk = Po^fc , where the polytropic variable 

^.^l+ p.^"^ ^ ^(00-0.) . (4) 

r(l + drel) Pq 

Here 6rei is a nonthermal component of pressure, assumed to be everywhere proportional to 
thermal pressure; the total is thus Ptot = (1 + 5rei)P- The pressure Pq and density po at the 
cluster center are found using two constraints: requiring conservation of energy, and matching 
the external surface pressure. For a given Pq and po, the final radius rj oi the gas initially 
inside r^r is found by summing outwards from the cluster center until the initial mass Mg is 
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enclosed. Note this means that gas may expand or contract, changing the gas fraction inside 
Tjjir, and also doing mechanical work. Assuming surface pressure changes little with radius, 
the change in energy is proportional to the change in volume: AEp = (47r/3)(rj^,y — r^)Ps- 



The constraint equation for conservation of energy is thus (see iPaper it ): 



Ef=Yl \po0f-'(pk + ^{l + 26rei)Po0f-' \l^ = Eg + AEp + eD\En\ + eFMpc^ . (5) 



rk<rf 



Here there are two additional inputs to the energy budget. The penultimate term represents 
energy input to the gas from dynamical processes, assumed to be proportional to E^', the 
final term represents feedback energy from collapsed objects, taken to be proportional to the 
mass of formed stars, Mp (different from M^, once gas recycling is taken into account — see 
Appendix lAl). These additional energy inputs will be described in Sec. [3l Finally, requiring 
the final surface pressure to match the external pressure yields the second constraint 



rf r 



:i + 6rel)Kf'Y.P^(^k-' =Ps , (6) 



k=l 



summing over the Nrf bins in the radial range rj < r < rf + I. (This is a change from 



Paper III , which used a wider zone, thereby giving a lower measurement of surface pressure). 



2.2. Creating the model cluster catalog 



The gas prescription was ap plied to c lust ers selected from a large-scale N-body simula- 
tion. We will refer the reader to lPaper 111 and ISehgal et al.l (120071 ) for further details, as the 
numerical methods used to carry out this simulation are the same as described in those pa- 
pers. The cosmological parameters were chosen to be {Qb, ^m, ^a, h, n.^^an) = (0.044, 0.264 



0.736, 0.71, 0.96, 0.80), consistent with the WMAP 5-year resuhs flKomatsu et al.l 120081 ). 
The simulation volume is 1000/i~^Mpc on a side and contains N=1024^ particles, making 
the particle mass 6.8 x lO^^/i^^M©; the particle spline softening length is 16.3/i~^kpc. At 
each simulation step particles were saved in one octant of a thin spherical shell centered on 
the origin, thus building up the matter distribution in a past light cone covering one eighth 
of the sky. All halos with a friends-of-friends mass of at least 2 x lO^^h'^MQ (293 particles) 
were selected from this light cone, out to redshift z = 0.2. Masses and radii at the virial 
and other overdensities are computed, in particular M500 and rsoo at a mean overdensity of 
500 times critical. The amount of star formation and iron production is determined based 
on the cluster redshift, as described in Appendix |Al The gas prescription was then run on 
each halo; there are over 15,000 halos in this catalog with M500 > 10^^h~^MQ. 
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To test that resolution is sufficient, a smaller sample of halos was selected from a 
higher resolution simulation. This has the same cosmological parameters, but the volume 
is 320/;,~^Mpc on a side, and the softening length is 3.2/i~^kpc; thus mass resolution is in- 
creased by a factor of 30.5, and spatial resolution by a factor of 5. No significant difference 
was seen when comparing results for lower mass halos (M500 ~ 1 — 3 x 10^^ h~^MQ) at the two 
resolutions; halos more massive than this are difficult to deal with at the higher resolution 
because of the large grid required. 

Once the density, temperature, and metallicity in each cell are kn own, the X-ray lu - 
minosity in the 0.5-2 keV band is computed w ith the M E KAL model (iMewe et al.lll985l ). 
For the cluster temperature, the mixT code of I VikhlininI (120061 ) was used to compute the 
spectroscopic temperature Tx in the projected annulus [0.15 — Ij-Rsoo! this is the definition 
of cluster temperature used throughout this paper. 

There are a number of parameters which must be specified: f^{Mmr),€F,eD, and 6rei- 
We will calibrate the model by comparing to X-ray data, in particular by matching the gas 
fraction as a function of Tx- Two subsamples from are selected from the halo catalog to 
match the properties of the two data sets we use for the calibrati on (this data is shown 
as a function of M500 in Figs. [T]and[2]). The ffist data set is from IVikhlinin et al.l (120061 ) . 
which consists mostly of higher Tx clusters; so, for comparison, we select all halos with 
-^500 > 3 X IO^^Mq in the z < 0.2 light cone. There are 76 halos in this subsample, with 
temperatures above 4 keV for most parameter choices. For a given model, we compute a 
Xv = '^[fg{model) — fg{obs)]'^/a{obs), summing for each observed cluster over all model 
points having a Tx within 4-cr of the observed T v . The second observational set is of low 
redshift groups (Tx ~ 1 — 2keV) from ISun et al.l (120091 ): so, for the second subsample, we 
select all model halos with z < 0.0485 and M500 > 3 x IO^^Mq. This subsample also contains 
76 halos, most with Tx < 3keV. A Xs is computed from this sample in an analogous manner 
as above. When searching for a best ffi model we will minimize = Xv + O-^xli the factor 
of 0.5 ensures that the two terms contribute roughly equally near the minimum value. 



3. The Role of Feedback and Star Formation 

It is instructive to ffist consider the case of no star formation or energy input, i.e. 
= /* = = 0; we will refer to this as the "Zero" model. The resulting temp e rature 



and Lx are shown as a function of mass as dotted lines in Fig. |H IVikhlinin et al.l (120091 ) 
found the M — Tx relation is well fit by a power law, E{z)M^qo = Mo(A;T/5keV)", with 
Mo = (3.02 ±0.11) X lO^^/i-^M© and a = 1.53. Fig. [1] shows the fractional difference {Tx - 
Tx{obs))/Tx{obs), where Tx{obs) is the value obtained by plugging M500 into the observed 
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relation. This figure shows that at a given M500, the predicted temperature from this simple 
model is lower than observed, by over 10%. At the same time, it shows that the predicted 
gas fraction remains near the cosmic mean, and is independent of mass (quite different from 
what is observed). This resuhs in a Yx = kTxMgir-^Qo) which is too high (by ~ 20% at higher 

r In 

masses) and an M^qq — Yx relation with a slope near the self-similar value Yx oc M^^^ , slightly 
shallower than observed. Thus we see the simplest set of assumptions yields an incorrect gas 
energy. There are three main processes which can alter this energy: transfer of energy from 
the DM through gravitational interactions, star formation, and feedback from astrophysical 
sources. In this section we will consider each of these in turn. 



So far (as in lPaper 111 ), it has been assumed that the gas originally had the same specific 
energy as the DM. However, this ignores an important process during cluster formation: when 
mergers occur, energy is transferred from the DM to the gas. In fact, the DM can lose up to 
10% of its energy in an equal mass collision, an d even a 10:1 merger can reduce the DM energy 
by a couple of percen t (IMcCarthy et al.ll2007l ). Simulations show such mergers are common: 
Cohn fc White! (120051) calculate that a ty pical cluster will experience four 5:1 accretion events 
since z ~ 2, and iFakhouri &: Mai (l2008l ) find that the mean rate for 10:1 or larger mergers 
in clusters is roughly two per halo per unit redshift. The magnitude of this effect can 
be estimated from hydrodynamical simulations, where it has long been found t h at baryon 



fraction in clusters is lower than the universal value 


(e.g. 


Cen & Ostriker 


1994; 


Eke et al. 


1998; 


Ettori et al. 


2006; 


Gottlober & Yeues 


2007; 


Stanek et al. 


2009). For 


simulations with 

average baryon 


no star formation or radiative processes. 


Grain et al.l ( 


200' 


'') find that the 



fraction within r2oo is 90% of the cosmic mean at 2=0. This result is independent of cluster 
mass, and at z=l the typical fraction is only 2 or 3% higher. Without star formation 
and energy feedback, the method of Sec. 12.11 gives a mean fg{< r2oo) ~ /c for halos with 
M200 > IO^^/i-^Mq (at low redshift, z < 0.2). Instead assuming that the energy transferred 
to the gas is 5% of the DM energy, i.e. us ing en = 0.0 5 in E qn. El this fraction is reduced 



to 0.92/c. As we are using a lower cxg than ICrain et al.l (120071 ) . a slightly higher gas fraction 
at z = seems appropriate. Thus we adopt en = 5% for the rest of the paper. 

Fig. [1], and also Fig. [2], display the change in fg inside v^qq after including this added 
energy, as short-dashed lines. While lower than before, it is still higher than observed, and 
still roughly independent of mass. The temperature at a given mass increases, but only 
slightly, so the M — Yx relation is still too high and shallow. Including a term proportional 
to Ed in this fashion does not change the self-similar nature of the model. 

Star formation will obviously decrease the gas fraction, but it has a second effect as well: 
since it is the most bound, lowest entropy gas that is turned int o stars, the remaining gas will 
have a higher specific energy than before (jVoit fc Bryanll200ll ). The simplest assumption is 
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to make the efficiency of star formation independent of cluster mass; Fig. [2] sliows, as dot- 
daslied lines, a = 0.05 model after also setting f^=9% for all halo masses. The predicted 
gas fraction drops further — for massive halos, near the observed value. Note the change in 
fg is more than 9%; the higher specific energy of the remaining gas is a significant effect, 
causing that gas to expand and further reducing the gas fraction in the inner parts of the 
halo. Temperatures also increase, so at higher masses the predicted Yx values are close to 
those observed. However, as the method up to this point has been scale-free, the gas fraction 
still does not depend on mass; thus, for less massive clusters, the predicted fg and Yx are 
still too large, since the slope is still the self-similar value. 

There are two ways to change the scale-free nature of the models: vary star formation 
efficiency with halo mass, or include feedb ack energy prop ortional to stellar mass. There 



is observational support for both options. iLin et al.l (120031 ). with a sample of 27 clusters 



measured the fraction of stellar mass to be larger for lower mass clusters, following M*(< 
^50o)/M5oo = 0.0164(M5oo/3 x IO^^Mq)""* with a, = 0.26 ± 0.09. To apply this to our 
model, we will assume this relation also holds inside the virial radius, i.e. M^/Md = M*(< 
'"50o)/j^50o at 2 = 0. A model varying with mass in this manner is shown in Fig. [1] as 
dot-dashed lines. As one would expect, this decreases the gas fraction more in lower mass 
halos; it also changes the slope of the M — Tx relation slightly. 

The slopes of the M — Tx and M — Yx relations wi ll clearly depend on how strongly 



varies with halo mass. Recently I Gonzalez et al.l (120071 ) found a much steeper relation in 



a sample of 23 clusters: M^,(< r5oo)/M5oo oc Mg^Q^^. A model where follows this relation 
is shown as long-dashed lines in Fig. [H The effect is quite dramatic in low mass clusters 
because so much gas is converted into stars (~ 50% or more), and the resulting T x and gas 



fractions tend to overshoot the observed trends at low masses. If we start with the Lin et al. 



( 120031 ) and vary the slope, searching for a best-fit model as described in Sec. 12.21 we find an 
exponent of a^^ = 0.49 yields a cluster sample with M — Tx and M — Yx relations quite close 
to those observed; this model is shown as solid lines in Fig. [TJ We will refer to this choice of 
parameters (e^ = 0.05, oc M^qq^^, ep = 0, 5rei = 0) as our "SF-only" model. 

The second effect which can break the scale-free nature of the models is additional 
energy injection, e.g. from SNe or AGN. We will assume such feedback is proportional to 
the number of stars formed, Mp (see Appendix |A]) . As an estimate of the amount of energy 
under consideration, for /*=9% and feedback energy lO'^Mpc^, the energy added amounts 
to 0.8 keV per particle initially inside r^r- Core collapse SNe could only add 0.4 keV/particle 
(assuming the model in Appendix |X] and 10^^ erg per event), so in this case there must be 
some additional source such as AGN energy from accreting black holes. If the mass in black 
holes is 10~^M*, and 2% of the rest mass energy of the accreted material is converted into 
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thermal energy (lAllen et al.ll2006l ). then an additional 1.0 keV/particle could be added to 
the gas. 

Fig. [2] shows the effect of adding feedback, starting with er, = 0.05 and = 9%. The 
long-dashed lines are for an extreme model with ep = 2.0 x 10^^. For the most massive halos, 
the feedback energy is not significant compared to the binding energy of the cluster; hence 
there is little change evident. For less massive clusters, the added energy is enough to cause 
the gas to expand, lowering the gas fraction considerably. The M — Tx relation becomes 
less steep as well, as lower mass halos become hotter; the M — Yx relation is thus less steep 
also, though in this example the change is so large as to be incompatible with observations. 

Varying the amount of feedback, the best-fit search yields a more moderate level of 
1.2 X IQ-^Mpc^, shown as solid lines in Fig. d We will call this {eo = 0.05, /, = 0.09, ep = 
1.2 X 10~^, 6rei = 0) the "Feedback" model. Like the SF-only model, this reproduces quite 
well the observed M — Tx and M — Yx relations. The efficiency 1.2 x 10~^ is smaller than 



the equivalent number used in lPaper 111 for several reasons. Part of the difference is merely 
bookkeeping: iPaper 111 did not consider mass loss from stars, and expressed feedback as 
being pro portional to M*, while here we instead use the larger mass Mp. More substantially. 
Paper III did not include the \Ed\ term, instead effectively considering it as part of the 
feedbac k. Finally , here the gas surface pressure is calculated over a smaller radial buffer zone 
than in iPaper III , which leads to a slightly lower gas fraction initially. 



Figs. [T] and [2] show that there is a degeneracy between steepening the star formation 
efficiency as a function of halo mass and increasing the amount of feedback per unit stellar 
mass. This is shown schematically in Fig. [31 We have found two models which both reproduce 
the X-ray data reasonably well: the SF-only m odel, which varie s s tar formation efficien cy as 



/* oc MgQQ (s teeper than the values found by lLin et al.l (120031 ) or lGiodini et al.l (120091 ) . but 



less steep than iGonzalez et al.l (120071 )). and the Feedback model, which instead has constant 
/* = 9% and adds feedback 1.2 x 10~^Mpc'^. A third option is of course some combination of 
these two. Since a stellar mass fraction independent of mass does appe ar to be contrad icted 
by observations, we will adopt the star formation efficiency as given by iLin et al.l (120031 ) and 
vary to find the best-fit model. This results in a lower level of feedback, ep = 4 x 10^^. 
We will call this our "Standard" model (eu = 0.05,/* oc Mgo^-^^ei. = 4 x 10'^, 6rei = 0), 
which we will use for the rest of the paper. Fig. |4] shows the M — Tx and M — Yx relations 
for this third model. This choice of paramete rs yields a popu l ation with thermal energy a s 
a function of mass quite close to that seen by IVikhlinin et al.l (120091 ) and ISun et al.l (120091 ) . 
However, keep in mind that there will be a range of models around this one which will fit 
the X-ray data equally well- it would be difficult to distinguish between these three models 
based solely on the gas profiles. The shaded region in Fig. [31 gives a rough idea of the 
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acceptable range of and ep- 



There is one further physical effect to be cons idered, namely nonthe rmal pressure sup- 



port, e.g. from cosmic rays or magnetic fields (e.g. iPfrommer et al.l 120071 ). We will conclude 
this section by considering how a moderate level of nonthermal support affects our model. 
If in fact 6rei = 0.10, then the M500 derived from X-rays are too low, because they are de- 
rived by assuming hydrostatic equilibrium. Fig. [5] shows our Standard model, after adjusting 
masses in the observational points by 10% to account for this; in this situation the model 
now has too high Tx and fg at a given M^qq. Adding 6rei = 0.10 to this model means that 
the gas requires less thermal pressure to maintain hydrostatic equilibrium, so the solution 
to Eqns. [5] and [6] yields gas which is cooler and correspondingly denser. This is shown as 
dashed lines in Fig. \5\ the model is now too cold, and the gas fraction even higher. Com- 
paring to Fig. [H it can be seen that increasing 6rei has the opposite effects of increasing e£). 
If we take the 6rei = 0.10 model and allow the dynamical energy input to vary, the best-fit 
result has to = 0.10, shown as solid lines in Fig. [51 This puts our predicted clusters back 
into agreement with the (adjusted) observations, at the same level as the Standard model. 
We will call this set of parameters (e^ = 0.10,/* oc Mg^o^*^, = 4 x 10~^,5rei = 0.10) the 
''Srei = 10%" model. 



Increasing e^j in this way assumes that non-radiative hydrodynamic simulations are 
missing some processes which transfer energy to gas during cluster formation. This certainly 
could be the case. For example, the dynamical friction experienced by the stars assernbling 
into a cD galaxy would add signific ant amounts of energy to the ICM ( lEl-Zant et al.ll2004j : 
Kimll2007l : 1Conrov fc Qstrikeilboosh . 



Radial Profiles 



Having fixed on a Standard model, in this section we will examine how well it matches 
observations at other radii than v^qq. 

Given that our DM halos have the standard density profile, the M2500 ~ Tx relation of 
our model should also agree with observations. Fitting all halos with kTx > IkeV, the best 
fit power law relation gives a normalization of M25ooE{z) = 1.25 x lO^^/i'^M© at kTx=5keV. 
This agrees well with obser ved X-ray samples, using masses derived fr om either weak lensing 
faoekstral [20071 ) or X-rays JArnaud et al.lbood : Ivikhlinin et al.lboOGh . 



Th e Standard model also has an X-ray luminosity in good agreement with the observa- 
tions of IVikhlinin et al.l (120091 ). unlike the Zero model which is a factor of 2-3 more luminous. 



This is hardly surprising, since we calibrated the model in part on this data. Above 10 M( 
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the predicted -^y oc M ^JJ , but below this the relation is much steeper, approaching M^(^^. 



Vikhlinin et al.l (120091 ) find the observed scatter at a given mass is ~ ±48%, while the model 
scatter is roughly half of this. There are many processes occurring in the cores of clus- 
ters which could significantly alter Lx not included in our model, which would increase the 
scatter and somewhat alter the mass dependence. 

A fundamental property of the gas is its entropy, usually defined in cluster studies as 
K = kT/rie^^, where Ue is the electron numj jer density. Based on Chandra observations 
of nearby relaxed clusters, iNagai et al.l (120071 ) determined how entropy varies with cluster 
temperature at several characteristic radii, including r5oo, tiooo, and r25oo- Their best-fit 
power-law approximations of these data are shown in Fig. [6] as dashed lines, which shows the 
entropy at these three radii as a function of Tx- The Zero model gives too low entropy values 
at all radii; true clusters have shallower entropy profiles than the scale-free model predicts. 
Once star formation and feedback are included, however, there is much improved agreement 
with obse rved clusters; th e main disagreement is at rsoo, where the model is roughly 10% 
too high. ISun et al.l (120091 ) also determined K — Tx relations at these radii, using Chandra 
archival data and including many lower temperature groups. For Tx > 2keV, the median 
value from the model cluste rs is within 10% o f their power- law fits at rsoo and riooo- At r25oo 
the model is lower than the ISun et al.l (120091 ) relation by roughly 15%, for all but the most 
massive clusters. Our model does not consider the history of cluster formation, nor the fact 
that the nonadiabatic processes which can alter the entropy will likely be more pronounced 
in the core. 



Afshordi et al.l (120071 ) assembled a heterogeneous catalog of 193 X-ray clusters and then 



examined WMAP data to derive a mean pressure profile using the SZ effect. This profile is 
shown as points with error bars in Fig. [7[ the pressure is normalized to the critical density 
times the global cluster Tx- For comparison, we computed the mean profile of the 284 
clusters from our Standard sample with kTx > 3keV. This is shown as solid line in Fig. [71 
the mean profile of the same halos without star formation or energy inp ut is shown as a dotte d 
line. To scale the radius, r2oo is estimated from kTx (as in Eqn. 1 of lAfshordi et al.l 120071 ). 
though using the radius measured directly from the DM distribution gives the same results. 
Outside of 0.5r2oo; changes to the gas energy show little effect on the mean profile. At smaller 
radii, star formation and feedback reduce the density more than they affect the temperature 
profile, thus lo w 'ering the pressure profile in the core. This gives better agreement with the 
Afshordi et al.l (120071 ) profile, though it is still higher. One source of this disagreement is 
that they assumed pT oc inside 0.25r2oo when constraining their profile, while our model 
is shallower than this; had they assumed a similarly shallow inner profile, this point would 
have been higher. Note also the redshift and temperature distributions of the model and 
observed samples are quite different. 
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Bonamente et al.l (120081 ) present a joint SZ and X-ray analysis of 38 clusters; here we 
will compare to their low-redshift sample of 19 clusters in the range 0.14 < z < 0.30. For 
our Stand ard sample with z < 0. 2 we calculated the integrated Compton y-parameter, Y 
(Eqn. 6 of iBonamente et al.ll2008l ). out to projected -R2500; th is is shown in Fig , p M ost of 
our sample is at lower gas masses than the clusters used by IBonamente et al.l (120081 ) , but 
where the two samples overlap the agreement is good. The best fit power-law found by 



Bonamente et al.l (120081 ). with slope 1.60, is shown as a dashed line in Fig. [HI A power-law 



fit to our model halos with Mg{< r25oo) > 8 x lO^^/i^^M© yields a slope of 1.61. Also shown 
as a dotted line is the Zero model. The effect of increasing the gas energy is to cause gas 
to expand and reduce the gas mass; however, this does not change the integrated pressure 
profile much, since the gas in the center must still hold up all the gas above it. This effect 
is larger for smaller mass halos, making the Y — Mg relation less steep. 

In Sec. [3] we found several different models which, given the degeneracy among input 
parameters, gave similar results. Since one goal of this is to provide a means of predicting 
the SZ signal from clusters, it is useful to compare these models. Fig. [9] shows the fractional 
difference in the Y value from the Standard model, as a function of Mg{< r25oo)5 for the SF- 
only. Feedback, and Srei = 10% models. The median SZ signal is within 10% of the Standard 
model in each case, and they are all roughly a factor of 2 higher than the Zero model, which 
does not include star formation and energy injection. This highlights the importance of 
normalizing to X-ray data when modeling the SZ signal. 



5. Summary and Discussion 

In this paper we have described an improved method for calculating the gas distribution 
in a DM potential well, constraining the model to match the observed X-ray gas fractions; 
then we made further comparisons with X-ray and SZ observations in order to check, where 
possible, our conclusions. Extracting the greatest return from upcoming cluster surveys will 
require understanding the state of the intracluster gas. It is encouraging that the simple 
model presented here can so well reproduce observed profiles and trends with cluster mass. 
With this model the following main points have demonstrated: 

• Simulations show that during mergers the energy of the DM component will be tapped to 
provide shock heating to the gas. Transferring ~ 5% of the binding energy of the DM to 
the gas results in gas fractions at large radii which are in good agreement with non-radiative 
hydrodynamical simulations. 

• Breaking the scale-free nature of the model to bring it into accord with observed scalings 
requires either varying the efficiency of star formation with cluster mass or increasing the 
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specific energy of tlie gas by an amount tliat is not proportional to tlie cluster binding energy 
(e.g. by being proportional to the stellar mass), or both. 

• Without including feedback energy, varying the efficiency of star formation as oc M, 



-0.49 



reprod uces observed M — Tx and M — Yx relations. This is steep er than found b v iLin et al. 



( I2OO3I ) and iGiodini et al.l (120091 ). but less steep than found by iGonzalez et al.l (120071 ). A 



relation as steep as oc M^qq^^, as seen by the latter, would leave too small gas fractions 
(and that gas at too high temperatures) in less massive clusters; however, the difference is 
less than two sigma. Also, keep in mind that what is measured is the stellar fraction inside 
rsoo, whereas /* in our model is the fraction inside Vmr; the two need not be identical. We are 
assuming that stars have the same radial distribution as D M, but if stars a r e inst ead more 
centrally located, this would change our predicted relation. iGonzalez et al.l (120071 ) selected 
for clusters with a dominant central galaxy, and for the least massive clusters r5oo is in the 
outskirts of this galaxy. Because a central dominant galaxy takes up a larger fraction of the 
stellar mass in smaller clusters, one could argue that the stellar fraction inside rsoo changes 
more rapidly with cluster mass than does the fraction inside Vmr, i.e. M^,(< rsoo) should vary 
with halo mass at a larger rate than 

• A model with M^/M^qq independent of cluster mass but including feedback energy pro- 
portional to the mass of formed stars can also reproduce the observed M — Tx and M — Yx- 
However, this requires adding roughly one keV per particle initially inside Vmr. This might 
be possible if the output from accreting black holes is converted to thermal energy with high 
enough efficiency (~ 0.02). 

• The most likely p ossibility is that both star formation and feedback are at work. Assuming 



oc M~°-^^ (as in iLin et al.l l2003l ) and including a lower feedback level of 4 x 10~^ Mpc^ 
also yields the proper scalings. Assuming the mass in black holes is 10~^M^,, this energy is 
6.5 X 10^'^ of their rest mass (or 0.33 keV per particle; note this does not take into account 
the feedback supplied by supernovae). This model was compared to further observations in 
Sec. m but the other two possibilities discussed here behave in quite a similar fashion. 
• Including nonthermal pressure results in a cooler, denser gas distribution. To bring such a 
model back into agreement with X-ray observations requires increased energy input. Once 
this is done, the pressure or entropy for a given Tx is quite similar to models with 6rei = 0, 
although of course the M — Tx relation is different. 

A couple caveats are in order. We start with halos drawn from a DM-only simula- 
tion; adding baryon s to such a simulation may change the properties of the halo population 
( IStanek et al.l 120091 ). Also, we have neglected many pertinent physical processes, such as 
cooling in the core (note our definition of cluster temperature excludes the core) . Because of 
this, the model halo population likely has too low a scatter at a given mass. This is exacer- 
bated by the tight relation we assume between halo mass and star formation plus feedback. 
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For example, if we allowed for a variation in at a given halo mass, clusters with higher 
would have higher temperatures and lower gas fractions, and thus lower X-ray luminosities. 
This would explain why c atalogs of optically selected clusters tend to be underluminous in 
X-rays at the same mass (iRykoff et al.ll2008l ). 



For most halos in our standard model, the gas initially assumed to be inside r^r has 
expanded outwards. For hotter {kT > 3keV) clusters only, the mean final radius Vf is 
1.23 ±0.13 (one sigma); for all clusters the mean Vf is 1.47±0.30. A result of this is that for 
the most massive clusters the baryon fraction inside is roughly 90% of the cosmic mean, 
fc = Qb/^m- Fig. [To] shows the the baryon fraction inside rsoo as a function of cluster mass 
for differing input parameters. The models shown have similar gas fractions (by design), but 
varying stellar mass fractions. This results in differing total baryon fractions at lower masses 
(though, as pointed out above, this may change if stars are more centrally concentrated than 
DM). In every case, however, the baryon fraction inside rsoo tops out near 0.8 fc at the highest 
cluster masses. Our sample contains four halos with virial masses above lO^^/i-^M©; their 
gas fractions within r2oo vary from 0.71/c to 0.82/c, with a mean of 0. 78/^. This agr e es we ll 
with the detailed mul tiwavelength analysis of four massive clusters by lUmetsu et al.l (l2008l ). 
Gonzalez et al.l (120071 ) also found b aryon fraction s inside rsoo near 0.8/c at high masses. As 
this work was nearing completion, iGiodini et al.l (120091 ) published a new determination of 
the baryon fraction (not including intra-cluster light) as a function of cluster mass. Our 
Standard model is in good agreement with their result, shown as a dotted line in Fig. [TOl 

The "missing" gas in this model is distributed in the outskirts of clusters. There it 
would be d i fficult to observe directly, being at low temperatures and densities. However, 
Prokhorovl (120081 ) has predicted in this case there would be emission in the extreme ul- 
traviolet from these outer shells. Alternatively, as this gas is associated w ith clusters, its 



kinetic SZ signature should be cor related with the cluster velocity field (IHo et al.l 12009 
Hernandez-Monteagudo fc Hdl2009l ). 



Computer simulations and analysis were supported by the National Science Foundation 
through TeraGrid resources provided by Pittsburgh Supercomputing Center and the National 
Center for Supercomputing Applications under grant AST070015; computations were also 
performed at the TIGRESS high performance computer center at Princeton University, which 
is jointly supported by the Princeton Institute for Computational Science and Engineering 
and the Princeton University Office of Information Technology. PB was partially supported 
by NSF Grant 0707731. 
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A. Star formation, gas recycling, and metal formation 



We will follow the "fossil" model of iNagamine et al.l (120061 ) . which matches local stellar 
populations; this represents a bulge and a disk component with delayed exponentials. Let 
Mp{Tj,t) denote the mass in stars of type j formed by time t, and Fj = Mp{Tj, tu) the total 
mass of such stars ever formed by a time tu corresponding to redshift zero, i.e. today. (At 



this point no allowance is made for mass lost from stars. ^ 
stellar mass we will adopt 



For the time evolution of this 



1+1 



This choice means the star formation rate is given by a delayed exponential 



MF{Tj,t) 



X*iTj,tH) 



X*{^j,t) 



1 t 



(Al) 



(A2) 



Nagamine et al.l (l2006l ) found a combination of a bulge population with decay time Tb=1.5 
Gyr and a disk population with rrf=4.5 Gyr is most consistent with local stellar populations 
(hence the term fossil). Since stellar populations in clusters are older than average, we will 
assume 90% of the total stars formed by ^ = are from the bulge component, and the 
remaining 10% from the younger disk component. The resulting star formation history is 
shown in Fig.[TTl there is little new star formation for z <1. This figure assumes M*/Mg=0.1 
at z=0; note that if increases then the amount of gas inside Vyij. will decrease, changing 
the final value of this ratio. 

The supernova rate also depends upon the delay time distribution. For a burst of star 
formation at t=0, suppose the fraction of SNe of type i occurring by time t is given by 
/5jv(Ti,t) = 1 — e~*/^\ Then the number of SNe per unit stel l ar mass per unit tina e will 
follow the commonly used distribution (e.g. IStrolger et al.ll2004j : iMannucci et al.ll2006l ) 



$(ri,t) = UifsNiTi] 



-tin 



(A3) 



where z/j is the total number of SNe of type i per unit mass which will ever occur. The SN 
rate S", i.e. the number of SNe per unit time, is then given by 

ft 



S{T„Tj,t) 



ds X*{'Tj,tH) 



fsN{Ti,t) + X*{Tj,t) 



Ti 



(^ 









(A4) 

where /ijj = Tj/(rj — Tj); we will be assuming Ti ^ Tj and Tj > throughout. Finally, the 
cumulative number of SNe of type i which have occurred by time t is 

dS{Ti,Tj,s) ^ 



S{Ti,Tj,t) 



ds 



-ds 



X*{rj,tH) 



J^iX{Ti,Tj,t) 



(A5) 
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fsN{ri,t) + 



(e-*/^^ - 1) 



(A6) 



Thus it only remains to specify Vi and for each type of SN. 



For core collapse SNe, the delay time is short; we will adopt tcc = 0.01 Gyr. We will 
employ a simple IMF with a slope o f -0.5 for stellar masses in t he range 0.1 — O.SMq and 
Salpeter slope -1.35 for 0.5 — IOOMq (IBaldry fc Glazebrookll2003l ). Assuming all stars in the 
range 8-50 M© end their lives as SNCC, then vcc ~ 9.14 x 10~^. 

The picture for SNIa is more complicated. There is growing evidence for a prompt 



SNIa component with a short delay time af t er star formation (e . g. iDella Valle et al. 



Mannucci et al.ll2005l . l2006l : iNeill et al.ll2006l : ISuUivan et al.ll2006l : lAubourg et al 



2005 



2007, and 



references therein). On the other hand, much observational evidence points towards much 



longer delay times of 2-4 Gyr (e.g. [Gallagher et al.ll2005l : iMannucci et al. 



2006 



2005 



Sullivan et al. 



Forster fc Schawinski 2008 . and references therein). Scannapieco fc Bildsten ( 20051 ) 



proposed that two evolutionary channels exist, a prompt population that would be propor- 
tional to the recent star formation rate, and a t ardy one with a much lo nger delay time that 
would be proportional to the total stellar mass. IMannucci et al.l (120061 ) found that the data 
are best matched by equal amounts of SNIa in two populations, one exploding of order 10^ 
yr after star formation and the other following an exponential distribution with a decay time 
of about 3 Gyr. We will follow this model for SNIa, setting a prompt delay time Tp = 0.05 
Gyr for one half, and a tardy Tj = 3 Gyr for the other half. 

Assuming SNIa result from 3-8 Mq stars, then Ui ^ 0.02754?7j, where rji is the fract ion of 
the st ars in this mass range that go supernova. This fraction is not well constrained; iMaoz 
( I2OO8I ) examines a number of different estimates of this fraction and finds that they vary 
from 2% to 40%, but also that r] ~ 15% is consistent with all the estimates he considered. 
Thus we will set rjp = rjt = 0.07, or Vp = ^1 = 0.00193, such that there are equal numbers of 
prompt and tardy SNIa. 

In a {h.VL^.K) = (0.7,0.3,0.7) cosmology, this model gives predicted S N rates per unit 
mass at z = of 0.07 SNuM for both SNCC and SNIa. For comparison, IMannucci et al. 
koO^ find 0.090l|J:S^^ and 0.070+[J;[J}^ SNuM (la statistical errors) for SNCC and SNIa 



respectively. The predicted rates increase with redshi ft; the predicted SN Ia rate is well 
within the la (statistical) errors of the rate r neasured bvlsharon et aL (2007) in the redshift 
range 0.06-0.19, and of the rate measured by I Graham et al.l (120081 ) in the range 0.2-1. From 
the supernova rate one can predict metallicity evolution. Using iron yields per event of 
O.O77M0 for SNCC and O.749M0 for SNIa, the resulting history of iron production is also 
shown in Fig. [TTJ Even though there is little star formation at low redshift, iron production 
continues because of tardy SNIa. Note the predicted metallicity will be proportional to the 
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amount of star formation; also note that we are assuming spatial homogeneity. Fig. [TT] also 
shows the predicted cumulative ratio of SNIa to SNCC; there are roughly 2.5 SNCC for each 
SNIa, and SNCC produce -20% of the iron. 

The rate of change in the actual stellar mass is the difference between the star formation 
rate and the rate at which mass is lost from stars. Codes exist for tracking the mas s loss 



from a stellar population, e.g. the PEGASE.2 code of iFioc fc Rocca-Volmerangd (119991 ). We 



will instead assume that all mass loss, whether from SNe or winds, can be modeled with 



Eqn. IA4[ Then the stellar mass M,, follows 

M,(r,-, t) = MpiTj, t)-Yl y^^^^'^ ^i' ^) ' (^'^) 

i 

where the mass yield Ui is the mean initial mass of stars of type i minus the mean remnant 
mass, if any. The delay times are adjusted to match the mass loss calculated by the PE- 
GASE.2 cod^ The cho,"ce= of delay times and yield s are summarized in Table [U We assume 
SNCC leave behind neutron stars of mass IAMq and SNIa leave no remnant. Stars larger 
than SOMq are all assumed to collapse to 2Mq black holes, on the same time scale as SNCC. 
Stars in the 5 — SM© range which do not go supernova instead remain as O.SSM© white 
dwarfs. The delay time for mass loss from these stars is set to 0.05 Gyr, slower than the 
prompt SNIa component; this choice is required to match the mass loss on short (< 1 Gyr) 
time scales. Stars in the 1.1 — SM© range also become white dwarfs, with mass loss on a 
much longer time scale of 6 Gyr; the lower limit of I.IM0 is treated as a free parameter, set 
to match the mass in white dwarfs at long time scales. Given these parameters, the evolution 
of the recycled gas fraction and the mass in remnants resulting from a single starburst is 
shown in Fig. [T21 also shown for the same IMF are the predictions of the PEGASE.2 code. 
Our relatively simple scheme follows the more detailed evolutionary code reasonably well. 



REFERENCES 

Afshordi, N., Lin, Y.-T., Nagai, D., & Sanderson, A. J. R. 2007, MNRAS, 378, 293 
Afshordi, N., Lin, Y.-T., & Sanderson, A. J. R. 2005, ApJ, 629, 1 

Allen, S. W., Dunn, R. J. H., Fabian, A. C, Taylor, G. B., & Reynolds, C. S. 2006, MNRAS, 
372, 21 

Arnaud, M. & Evrard, A. E. 1999, MNRAS, 305, 631 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 

— . 2007, A&A, 474, L37 



^Available at [http : //www2 . iap . f r/users/f ioc/PEGASE . html 



-18- 



Ascasibar, Y. & Diego, J. M. 2008, MNRAS, 383, 369 

Ascasibar, Y., Yepes, G., Miiller, V., & Gottlober, S. 2003, MNRAS, 346, 731 

Aubourg, E., Tojeiro, R., Jimenez, R., Heavens, A. F., Strauss, M. A., & Spergel, D. N. 

2007, ArXiv e-prints, 707 

Babul, A., Balogh, M. L., Lewis, G. F., & Poole, G. B. 2002, MNRAS, 330, 329 
Baldry, I. K. & Glazebrook, K. 2003, ApJ, 593, 258 
Balogh, M. L., Babul, A., & Patton, D. R. 1999, MNRAS, 307, 463 
Bhattacharya, S., di Matteo, T., & Kosowsky, A. 2008, MNRAS, 389, 34 
Bode, P., Ostriker, J. P., Weller, J., & Shaw, L. 2007, ApJ, 663, 139 
Bohringcr, H. ct al. 2007, A&A, 469, 363 

Bonamentc, M., Joy, M., LaRoque, S. J., Carlstrom, J. E., Nagai, D., & Marrone, D. P. 

2008, ApJ, 675, 106 

Borgani, S., Diaferio, A., Dolag, K., & Schindler, S. 2008, Space Science Reviews, 134, 269 
Borgani, S., Dolag, K., Murante, G., Cheng, L.-M., Springel, V., Diaferio, A., Moscardini, 

L., Tormcn, G., Tornatore, L., & Tozzi, P. 2006, MNRAS, 367, 1641 
Bower, R. G., McCarthy, 1. G., & Benson, A. J. 2008, MNRAS, 390, 1399 
Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 

Burenin, R. A., Vikhlinin, A., Hornstrup, A., Ebeling, H., Quintana, H., & Mescheryakov, 

A. 2007, ApJS, 172, 561 
Cen, R. & Ostriker, J. P. 1994, ApJ, 429, 4 
Ciotti, L. & Pellegrini, S. 2008, MNRAS, 387, 902 
Cohn, J. D. & White, M. 2005, Astroparticle Physics, 24, 316 
Conroy, C. & Ostriker, J. P. 2008, ApJ, 681, 151 

Grain, R. A., Eke, V. R., Frenk, G. S., Jenkins, A., McCarthy, I. G., Navarro, J. F., & 

Pearce, F. R. 2007, MNRAS, 377, 41 
Dahle, H. 2007, ArXiv Astrophysics e-prints 

Dave, R., Oppenheimer, B. D., & Sivanandam, S. 2008, MNRAS, 391, 110 

Delia Valle, M., Panagia, N., Padovani, P., Cappellaro, E., Mannucci, F., & Turatto, M. 

2005, ApJ, 629, 750 
Dobbs, M. et al. 2006, New Astronomy Review, 50, 960 
Dos Santos, S. & Dorc, O. 2002, A&A, 383, 450 
Eisenhardt, P. R. M. et al. 2008, ApJ, 684, 905 
Eke, V. R., Navarro, J. F., & Frenk, C. S. 1998, ApJ, 503, 569 
El-Zant, A. A., Kim, W.-T., & Kamionkowski, M. 2004, MNRAS, 354, 169 
Ettori, S., Dolag, K., Borgani, S., & Murante, G. 2006, MNRAS, 365, 1021 
Fakhouri, O. & Ma, G.-P. 2008, MNRAS, 386, 577 
Finoguenov, A. et al. 2007, ApJS, 172, 182 

Fioc, M. & Rocca-Volmerange, B. 1999, ArXiv Astrophysics e-prints 
Forster, F. & Schawinski, K. 2008, ArXiv e-prints, 804 



-19- 



Gallagher, J. S., Garnavich, P. M., Berlind, P., Challis, P., Jha, S., & Kirshner, R. P. 2005, 

ApJ, 634, 210 
Giodini, S. et al. 2009, ArXiv e-prints 

Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 147 

Gottlober, S. & Yepes, G. 2007, ApJ, 664, 117 

Graham, M. L. et al. 2008, AJ, 135, 1343 

Hernandez-Monteagudo, C. & Ho, S. 2009, ArXiv e-prints 

Ho, S., Dcdco, S., & Spergel, D. 2009, ArXiv e-prints 

Hoekstra, H. 2007, MNRAS, 379, 317 

Kaiser, N. 1991, ApJ, 383, 104 

Kay, S. T., da Silva, A. C, Aghanim, N., Blanchard, A., Liddle, A. R., Puget, J.-L., Sadat, 

R., & Thomas, P. A. 2007, MNRAS, 377, 317 
Kim, W.-T. 2007, ApJ, 667, L5 
Koester, B. P. et al. 2007, ApJ, 660, 239 
Komatsu, E. et al. 2008, ArXiv e-prints 
Komatsu, E. & Scljak, U. 2001, MNRAS, 327, 1353 
Kosowsky, A. 2006, New Astronomy Review, 50, 969 
Lapi, A., Cavaliere, A., & Menci, N. 2005, ApJ, 619, 60 

Lin, K.-Y. et al. 2008, in Astronomical Society of the Pacific Conference Series, Vol. 399, 
Astronomical Society of the Pacific Conference Series, ed. T. Kodama, T. Yamada, 
& K. Aoki, 384- 

Lin, Y.-T., Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749 

Locwenstein, M. 2000, ApJ, 532, 17 

Lopes, P. A. A., de Carvalho, R. R., Kohl-Moreira, J. L., & Jones, C. 2009, MNRAS, 392, 
135 

Malte Schafer, B. & Bartelmann, M. 2007, MNRAS, 377, 253 

Mannucci, F., Delia Valle, M., & Panagia, N. 2006, MNRAS, 370, 773 

Mannucci, P., Delia Valle, M., Panagia, N., Cappellaro, E., Cresci, G., Maiolino, R., Pet- 

rosian. A., & Turatto, M. 2005, A&A, 433, 807 
Mannucci, F., Maoz, D., Sharon, K., Botticella, M. T., Delia Valle, M., Gal- Yam, A., & 

Panagia, N. 2008, MNRAS, 383, 1121 
Maoz, D. 2008, MNRAS, 384, 267 

McCarthy, L G., Bower, R. G., Balogh, M. L., Voit, G. M., Pearce, F. R., Theuns, T., Babul, 

A., Lacey, C. G., & Frenk, C. S. 2007, MNRAS, 376, 497 
Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197 
Moodley, K., Warne, R., Gohccr, N., & Trac, H. 2008, ArXiv c-prints 
Muanwong, O., Kay, S. T., & Thomas, P. A. 2006, ApJ, 649, 640 
Muchovej, S. et al. 2007, ApJ, 663, 708 

Nagai, D., Kravtsov, A. V., & Vikhhnin, A. 2007, ApJ, 668, 1 

Nagamine, K., Ostriker, J. P., Fukugita, M., & Cen, R. 2006, ApJ, 653, 881 



-20- 



Neill, J. D. et al. 2006, AJ, 132, 1126 

Ostrikcr, J. P., Bode, P., & Babul, A. 2005, ApJ, 634, 964 
Pacaud, F. et al. 2007, MNRAS, 382, 1289 

Pfrommer, C, Enfilin, T. A., Springel, V., Jubelgas, M., & Dolag, K. 2007, MNRAS, 378, 
385 

Pratt, G. W., Croston, J. H., Arnaud, M., & Boehringer, H. 2008, ArXiv e-prints 
Prokhorov, D. A. 2008, A&A, 492, 651 

Puchwcin, E., Sijacki, D., & Springel, V. 2008, ApJ, 687, L53 

Rykoff, E. S., McKay, T. A., Becker, M. R., Evrard, A., Johnston, D. E., Koester, B. P., 

Rozo, E., Sheldon, E. S., & Wechsler, R. H. 2008, ApJ, 675, 1106 
Scannapieco, E. & Bildsten, L. 2005, ApJ, 629, L85 
Sehgal, N., Trac, H., Huffenberger, K., & Bode, P. 2007, ApJ, 664, 149 
Sharon, K., Gal- Yam, A., Maoz, D., Fihppenko, A. V., & Guhathakurta, P. 2007, ApJ, 660, 
1165 

Shimizu, M., Kitayama, T., Sasaki, S., & Suto, Y. 2004, PASJ, 56, 1 
Short, C. J. & Thomas, P. A. 2008, ArXiv e-prints 

Sijacki, D., Springel, V., di Matteo, T., & Hcrnquist, L. 2007, MNRAS, 380, 877 

Solanes, J. M., Manrique, A., Gonzalez-Casado, G., & Salvador-Sole, E. 2005, ApJ, 628, 45 

Stanek, R., Rudd, D., & Evrard, A. E. 2009, MNRAS, L172+ 

Staniszewski, Z. et al. 2008, ArXiv e-prints 

Strolger, L.-G. et al. 2004, ApJ, 613, 200 

Sullivan, M. et al. 2006, ApJ, 648, 868 

Sun, M., Voit, G. M., Donahue, M., Jones, C., Forman, W., & Vikhlinin, A. 2009, ApJ, 693, 
1142 

Suto, Y., Sasaki, S., & Makino, N. 1998, ApJ, 509, 544 
Tozzi, P. & Norman, C. 2001, ApJ, 546, 63 
Umetsu, K. 2008, ArXiv e-prints 
Vikhlinin, A. 2006, ApJ, 640, 710 

Vikhlinin, A., Burenin, R. A., Ebeling, H., Forman, W. R., Hornstrup, A., Jones, C., 
Kravtsov, A. V., Murray, S. S., Nagai, D., Quintana, H., & Voevodkin, A. 2009, 
ApJ, 692, 1033 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C., Markcvitch, M., Murray, S. S., & Van 

Speybroeck, L. 2006, ApJ, 640, 691 
Voit, G. M. & Bryan, G. L. 2001, ApJ, 551, L139 

Voit, G. M., Bryan, G. L., Balogh, M. L., & Bower, R. G. 2002, ApJ, 576, 601 

Wu, K. K. S., Fabian, A. C, & Nulsen, P. E. J. 2000, MNRAS, 318, 889 

Yee, H. K. C., Gladders, M. D., Gilbank, D. G., Majumdar, S., Hoekstra, H., & Ellingson, 

E. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 379, Cosmic 

Frontiers, ed. N. Metcalfe & T. Shanks, 103- 



-21 - 



This preprint was prepared with the A AS IM^jX macros v5.2. 



Table 1. Decay times and yields 



mass 


mean 


remnant 


remnant 








n 


range 


mass 


type 


mass 


(fraction) 






(Gyr) 


50-100 


68.4 


BH 


2.0 


1 


5.11 X 10" 


-4 


0.01 


8-50 


16.0 


SNCC+NS 


1.4 


1 


9.14 X 10" 


-3 


0.01 


3-8 


4.58 


prompt lA 





0.07 


1.93 X 10" 


-3 


0.05 


3-8 


4.58 


tardy lA 





0.07 


1.93 X 10- 


-3 


3.01 


3-8 


4.58 


WD 


0.55 


0.92 


0.0253 




0.5 


1.1-3 


1.69 


WD 


0.55 


1 


0.1078 




6.0 



Note. — Masses in Mq. 
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Fig. 1. — Effects of vary ing star formation, Fractional difference from tlie observed mass- 



temperature relation of IVikhlinin et al.l (l2009l ). h) Stellar (thinner lines) and gas (thicker 
lines) fractions, normalized to the cosmic average. Shown are models with: no star formation 
and e/) = o r 0.05 (dotted a nd short-dashed lines, respectively); als o en = 0.05 and oc 



Lin et al.ll2003f ). oc M-°-^^ and /, oc M 



(as in 

(dot-dashed , solid, and long -dashed lines), 
circles) and ISun et al.l (120091, crosses) . 



-0.64 



,as m 



Gonzalez et a 



Data points are from IVikhlinin et al. 



2007) 



(12006 




Fig. 2. — Effects of varying feedback. As in FigH], showing models with: no star formation 
and ez) = or 0.05 (dotted and short-dashed lines, respectively); also eo = 0.05 plus fraction 
/* = 9% of gas converted to stars, with feedback energies of 0, 1.2, and 2.0 x l{)~^Mpc^ 
(dot-long dashed, solid, and long-dashed lines). 
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Fig. 3. — Schematic depiction of the degeneracy between varying star formation efficiency 
/* oc M^Q* and feedback energy epMpc^. Symbols show three models based on different 
assumed /*, as marked. The shaded region denotes a constant Ax^ from the Standard 
model, corresponding to a 99% confidence interval for two degrees of freedom (this is only 
approximate, as we have not minimized over other input parameters). 



-24- 




Fig. 4. — Mass-temperat ure and mass-Yx relations for the Standard model with en = 0.05, 
/* as in iLin et al.l (120031 ) . and ei? = 4 x 10~^. Sohd lines are the median values, and hash 
marks denote the region enclosing 68% of the halos. At the high mass end individual halos 
are shown as x's. Dotted lines give the median for the Zero model (no st ar formation or 



feedback). Straight lines are the observed bes t-fit power laws: dashed from I Vikhlinin et al. 



( I2OO9I ) ar id dot-dashed from ISun et al.l (120091 ). Lower panels show the fractional difference 
from the I Vikhlinin et al.l (120091 ) relations. 




Fig. 5. — Effects of nonthermal pressure. As in Fig{T], but with observational data adjusted 
to assume a 10% underestimate of cluster mass. Shown are: the Standard model (dotted 
lines), after including 5rei = 0.10 (dashed lines), and also doubling eo (solid lines). 
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Fig. 6. — Entropy as a function of tempera ture, at r^sn n , ^innn , and t^qq. Straight dashed 
hues are the observed best-fit power laws of iNagai et aLl (120071 ). Sohd hues are the median 
values for the Standard model, dotted lines for the Zero model. For the former, clusters at 
high Tx are sho wn as individual p oints. Lower panels shows the fractional difference of the 
model from the iNagai et al.l (120071 ) relations. 
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Fig. 7. — The SZ profile as a function of radius. Solid line is the mean profile for all 
/cTj, > 3keV clusters; dashed lines show two standard deviations. The dotted shows the 
mean profile for the Zero model. Po ints with error bars are the mean profile found from 
WMAP data bv kfshordi et all toO^ . 




Fig. 8. — The SZ signal as a function of gas mass at overdensity 2500. The solid line is 
the median values, with hash marks denoting the region enclosing 68% of the halos; at 
the high mass end individual halos are shown as x's. The dotted line is the median for a 
model with no st a r form ation or feedback. The long-dashed line is the observed relation of 



Bonamente et al.l (l2008l ). 
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Fig. 9. — Fractional difference of the SZ signal Y from the Standard model. Shown are 
the Zero (dotted line), SF-only (short-dashed), Feedback (long-dashed), and 6rei = 10% 
(dot-dashed) models. 



V 

o 

o 

cd 

Sh 



0.8 



0.7 



0.6 



0.5 



0.4 





1 


1 1 1 


1 1 L 


































" / / 
/ / 
































/ / 
/ / 








_ / / 








/ / 
/ / 










1 




1 1 1 



1014 



101 



Fig. 10. — Baryon fraction inside rsoo as a function of mass. The lower set of lines represents 
gas fraction, and the upper set total baryon fraction. Shown are the Standard (solid lines), 
SF-only (short-dashed), and Feedback (long-dashed) mo dels. The dott e d line is the baryon 
fraction (not including intra-cluster light) determined by iGiodini et al.l (120091 ). 
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Fig. 11. — Star formation history as a function of redshift. Long dashed: mass in stars, as 
a fraction of all baryons. Solid: total mass of stars formed. Dotted: metallicity with respect 
to solar. Short dashed: the ratio of SNIa to SNCC. All quantities are measured inside rj, 
and assume the redshift zero star/gas ratio is M*/Mg=0.10. 




Fig. 12. — Fraction of gas recycled and locked up in remnants as a function of time for a 
single starburst at t = 0. From top to bottom, the lines are fraction of recycled gas (solid), 
in white dwarfs (dotted), and in black holes or neutron stars (dashed). For comparisoii , the 
same quantities as predicted by the PEGASE.2 code (IFioc fc Rocca-Volmerangelll999l ) are 
shown as circles, triangles, and squares, respectively. 



